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Translocation properties of ionic channels are investigated, on the basis of classical electrostatics, 
with an emphasis on asymptotic formulas for the potential and field associated with a point charge 
in the channel. Due to image charges in the membrane, we show that ions in an infinite length 
channel interact via a one-dimensional (fD) Coulomb potential. The corresponding electrostatic 
barrier E is characterized by a "geometric mean" screening E oc e^ j y/e-m^mR (R being the radius 
of the pore, and em ~ 2 and em ~ 80 the room temperature diele ctric co nstants of membrane 
and water, respectively). There exists a crossover length, xq oc i?^ ew/em ~ 6.3R, below which 
the ID potential governs the electrostatics and beyond which the three-dimensional (3D) Coulomb 
potential screened by the membrane takes over. Knowledge of this length enables us to discriminate 
between long channels, the length L of which satisfies: L S> 2xo, and short channels for which 
L <^ 2xo. The latter condition is satisfied by most realistic channels {e.g., gramicidin A where 
R « SA, L ~ 2.5nm and 2xo ~ 3.8nm) whose translocation energy is therefore controlled by the 
part of the self-energy, E, arising from the fD potential. On this basis, we derive an expression for 
E, with no fitting parameter, which applies to a generic nano-channel of length L and radius R. 
Our results are related to model-independent translocation properties of nano-scale ionic channels, 
they improve on previous, curve-fitting, formulas and agree to within 5% with estimates, resulting 
from numerical simulations, available in the literature on the subject. 



I. INTRODUCTION 



The simplest physical models of bilayer membranes 
turn out to be particularly useful to understand the basic 
electrostatic interplays between the membrane and the 
ionic flow. Consider the phospholipid bilayer as a thin 
slab of non-conducting material separating two aqueous 
solutions, and an aqueous ion channel as a cylindrical 
region connecting the two solutions, see Fig^ One ap- 
parent obstacle to ionic conduction is the electrostatic en- 
ergy barrier which has to be passed in order to move the 
ion through the channel. In 1969, Parsegiani found that, 
for an infinite channel, the electrostatic self-energy E, i.e. 
the energy required to move an ion from infinity to the 
center of the channel, is of the order of 16 ksT. A major 
contribution to this self-energy comes from the self-image 
of the charge due to the consequent difference between 
the dielectric constant of the water channel (e^, ~ 80, 
at room temperature) and that of the lipid membrane 
(ern ~ 2). Such a large barrier would prevent any ionic 
flow, contrary to the relatively high conductances ob- 
served in experiments^. Moreover, these are compatible 
with barriers not exceeding 6 fcsT. The most obvious 
limitation of previous models is that they consider chan- 
nels of infinite length (no screening by counter ions is 
assumed because typical channel radii considered here 
are of the order of a few A). Taking into account the 
finite length of the channel (gramicidin A has a charac- 
teristic length of 25 A), as well as a high dielectric shield 
of the channel, the authors of Refs. 3 and 4 have shown, 
with the help of numerical methods, a considerable de- 
crease of this barrier, reaching 6.7 ksT. A large number 
of recent studies also address this question, generally fo- 
cusing on specific models and taking into account of the 
channel nano-structure, using numerical simulations of 













L 






• 













X 



FIG. 1: Schematic view of a channel of length L, radius R 
with an ion sitting in its center, i.e. at a; = 0. The channel 
consists of water, of dielectric constant e^ ~ 80, and is sur- 
rounded by the membrane, of dielectric constant £,„ « 2. A 
typical channel we are interested in, e.g. gramicidin A, has a 
length L « 2.5nm and a radius R ~ SA. 



the Poisson-Boltzmann or Poisson-Nernst-Planck type, 
to improve on previous estimates^. 

In this contribution, the original, model-independent, 
problem of the electrostatics of point-like ions in an infi- 
nite channel is reconsidered, with an emphasis on asymp- 
totic expressions for the potential and field which have 
not been derived previously. Our quantitative results 
show the appearance of a length scale: xq cx R^ew/^m, 
which separates regions, in space, where the effect of 
the image-charges on the ion at the origin are cru- 
cially different. At distances x, along the channel, 
smaller than xq (x <^ Xq), the electrostatic potential, 
if, is one-dimensional (ID) and displays an unusual, i.e. 
algebraic in the dielectric constants, screening: ip oc 
{e"^ /y/eiue,nR){l—\x\/xo). At distances x, along the chan- 
nel, larger than xq {x ^ xq), the electrostatic potential is 
three-dimensional (3D) and everywhere, even in the wa- 
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ter channel, screened by the membrane: </? oc e^/em|a;|. 
In the region x Xq, the ions are anti-confined, i.e. sub- 
ject to a constant repulsive force. The crossover length, 
xq , is found to be of the order of, or larger than, the aver- 
age length of a typical ionic channel and our asymptotic 
expressions, in the ID Coulomb regime, thus capture the 
essence of the phenomena characteristic of realistic, i.e. 
nano-length, channels. On this physical basis, we will 
therefore derive a formula for S depending on the length 
L and radius i? of a typical pore. This formula for E, 
will be shown to agree, within 5% of accuracy, with the 
available numerical results in the field and therefore im- 
proves on previous, curve-fitting, expressions. This will 
lead us to an accurate knowledge, of the main exponen- 
tial dependence, of the conductance G, of finite-length 
channels: G oc exp[—'E,/kBT]. 

This article is organized as follows. We first consider 
the case of an infinite channel where calculations may be 
carried out straightforwardly. In Sec.^ the electrostatic 
potential, within the infinite channel, is derived in var- 
ious asymptotic regimes. The length xq is defined. In 
Sec. mil similar asymptotic formulas for the potential in 
the membrane are derived. These results are used, in 
Sec. IIVI to derive the components of the electric field. 
In Sec. the image self-energy of a charge in an infi- 
nite channel is derived and the part of this self-energy, 
arising from the ID Coulomb potential, is singled out. 
In Sec. IVII numerical estimates of the various quantities, 
derived in the previous sections, are given and compared 
with exact numerical simulations. In Sec. IVIII the case 
of finite size channels is considered and, on the basis of 
the results of previous sections, an expression for the self- 
energy of nano-channels is derived. A comparison is then 
made with the available results in the literature, based 
on numerical simulations. Finally, the conclusion is given 
in Sec. Ivml 



II. ELECTROSTATIC POTENTIAL INSIDE AN 
INFINITE CHANNEL 

Consider a point-like ion (charge e) in an infinite cylin- 
der of radius R and of dielectric constant Cw (of water), 
embedded in a media (the membrane) of dielectric con- 
stant Cm ^ iw, cf. Fig. ^ The electrostatic potential, 
ip{x,r±), of such an ion within the cylinder is given by: 



ip{x,r^) 



dk 



ip{k, r±) cos{kx), 



(1) 



where k is the wavenumber along the (x-)axis of the chan- 
nel and: 
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2e (e^ - ern)KoikR)KiikR)Ioikr^) 



(2) 



e„ e^Ko{kR)Ii{kR) + emKi{kR)h{kR) 

where J,^ and Ky are modified Bessel functions of the 
first and second kind, respectively. The result of Eqs. ^ 



and jSJ may be straightforwardly derived by solving the 
appropriate inhomogeneous Poisson equation. The first 
term in Eq. ^ corresponds to the usual (l/r in 3D) 
Coulomb potential of a charge in an infinite media of 
dielectric constant and the second term corresponds 
to the potential due to the image charges (repulsive, as 
it should be, for > Cm, because the charge is sitting in 
the high dielectric constant media) . A similar expression 
has been derived elsewhere^ and we shall therefore not 
prove it here. 

At distances smaller than the radius of the channel, 
kR ^ 1, we may use the asymptotic expressions of mod- 
ified Bessel functions^: 



K,{kR) 



exp(— fci?) 
y/2TTkR ' 



L{kR) 



exp{kR) 



(3) 



Substituting these expressions in the second term, the 
image part, of Eq. (0), leads to a term « e/e„i?, which, 
at X <C i?, is negligible compared to the first one: e/ew\x\. 
At distances, x <C i?, the images play no role, and the 
Coulomb potential is the usual 3D one, screened by the 
dielectric constant of water: 



e 



< R. 



(4) 



Because all phenomena we shall be interested in, in the 
following, are related to the effect of image charges, R will 
be our lowest distance cut-off. We will therefore explore 
the properties of the system at distances larger than the 
radius of the channel. 

We focus on the asymptotics of Eq. (|2Jl at distances, 
along the channel, larger than its radius, i.e. kR <^ 1. 
We will need, in the following, the series expansions of 
Bessel functions, which may be found in any standard 
table, e.g. Ref. 0- To the highest order of main interest 
to us, these series expansions read: 

KoikR)^\og{2/kR)-j, 

Ki(kR) « -(log(2/fci?) - i)kR/2, 

Io{kR)Kl + {kR/2f + 

Ii{kR) K kR/2 + (5) 

where 7 w 0.577 is Euler's constant. With the help of 
Eq. (|SJ), Eq. jSJ) may then be approximated by: 



(p{x,r±) 



ew\x\ 



+00 



dUKo{OK,iOcos{^x/R) 



-^{1- — ) 
^niR Jo TT 1 + {^xo/Rf 



(6) 



where the dimensionless ^ = kR has been introduced and 
a characteristic distance, xq, appears and reads: 



Xq = R 



2e„ 



Vlog(2a:o/i?) - 7- 



(7) 



For clarity, the solutions of Eq. {Tj) for the dimensionless 
ratio xq/R as a function of t^,l2trn are plotted on Fig.|21 
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FIG. 2: The crossover length xq, in units of the radius of the 
channel R, as a function of the ratio em/2em- 
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FIG. 4: Schematic view of an infinite channel, including the 
cylinder of length 2xo and radius so, cf. Eq. Q. 
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FIG. 3: The crossover length xo, in units of the radius of the 
channel i?, as a function of the ratio ew/2em in the vicinity of 
the realistic value e„,/2em = 20, for a water channel in a lipid 
membrane {tw ~ 80 and em ^ 2). 



Fig. 121 provides more accurate values for the crossover 
length xq in the vicinity of the reahstic value ew/2em = 
20, for a water channel in a lipid membrane {ew ~ 80 and 
em - 2). 

In the limit ^ <C -R/xq, i.e. x ^ Xq, Eq. 10 reduces 
to: 
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ip{x,rj_ 
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The integral is a standard one and, as a result 

(p{x,rj_) 



1 3^ I I I I I 



(9) 



which corresponds to the usual 3D Coulomb interaction 
screened by the membrane. 

In the limit ^ > R/xq, i.e. R ^ x ^ xq, Eq. © 
reads: 

g 

(y9(x,r_L) = i4:o(i?/xo)Jo(r±/a;o)exp(-|a;|/a;o), 

£mXQ 

i? < X < Xo, r_L < i?, (10) 

where corrections in x/xq, r±/xo as well as R/xq have 
been summed, to all orders, into the special functions 
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FIG. 5; The electrostatic potential at the center of an infinite 
channel {r± = 0) as a function of x/R, the dimensionless 
distance along the water channel {e^ ~ 80 and « 2). 
The tick at xo — 6.3R marks the crossover between the ID- 
potential, X <^ Xo, and the 3D potential, x xq. 



appearing in Eq. (|10|l . With the help of series ex- 
pansions of these special functions, to the lowest or- 
ders in the small quantities above: Joifj^/xo) ~ 1 and: 
Ko{rj_/xo) « -7-l-log(2xo/i?), Eq. ^ yields: 



^{x, r^) « ^ (log(2xo/i?) - 7) 1 - ^ 



i? < X < Xq, r± < R. 



(11) 



Eq. Hll() shows that the potential is ID-like, because of 
the linear dependence on the coordinate, x, along the 
channel. The electrostatic potential along the axis of the 
channel has been plotted on Fig. [S] from Eqs. (Q) and 



III. ELECTROSTATIC POTENTIAL OUTSIDE 
THE INFINITE CHANNEL 

In this section, we will still focus on the ion within the 
pore but determine the potential due to this ion outside 
the pore, i.e. for r± > R. This potential reads: 



2e 



kR e.^Ko{kR)h{kR) + e,„Xi(fci?)/o(fci?) 



(.12) 
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Being interested in distances, along the axis of the pore, 
which are larger than i?, i.e. kR <^ 1, Eq. H12I) reads, in 
this limit: 



ip{x, rj 



2e 



^m^O Jq 



+00 



Ko{^r±/xo) cos(ga:/a:o) 

TT 1+e 



(13) 



where the dimensionless variable ^ = kxo and xq is given 
by Eq. O. 

For X ^ Xq, we recover the potential of a 3D point-like 
charge: 

(p{x,r±) = , a;>a;o, V r_L > i?, (14) 



where r — a/x^ + and the screening is by the dielectric 
constant of the membrane. 

For a; ^ xo, the potential reads: 



, , 0/1 



[i/o(?'±/xo) - Yo{r±/xo)] 



(15) 



where iJo is a Struve function and Yq is a Neumann 
function. Using asymptotic and series representations 
of these special functions^, we find that: 



(p{x, r^) w — [log(2a:;o/r_L) - 7] , 
x ^ xo, R < r± xo, 



and that: 



ip{x,r^) = 



X <^ Xq, r± > Xq. 



(16) 



(17) 



Eq. (|16|l shows that, inside the shell of radius xq outside 
the pore, the potential is 2D-like. As for the ID potential 
of Eq. 1)11(1 , this implies that field lines are mostly located 
within the water channel because of its high dielectric 
constant. Beyond xq, such a potential crosses over to 
the 3D Coulomb potential of Eq. H17|l , similarly to what 
happens inside the channel, cf. Eq. 

The potential outside the pore of Eq. Hl()|l has also an 
x— dependence which we have neglected as a first approx- 
imation. Taking this dependence into account, the latter 
may formally be expressed as: 



ip{x,rj_) = — —Koir_i_/xo)ex.p{-\x\/xo), 
X Xo, R < r± <^ Xq. 



(18) 



The series expansions of the special functions appearing 
in Eq. (|18|) , give the dominant contributions in the small 
parameters x/xq and r±/xo, e.g. the exponential gives 
the linear dependence on |x| and the MacDonald function 
Kq gives the logarithmic dependence on r±. 



IV. ELECTRIC FIELD IN AN INFINITE 
CHANNEL 

The previous sections have shown the importance of 
the crossover length, xq, which separates regions, in 



space, characterized by different sectors of the electro- 
static potential controlled by the effect of image charges. 
We proceed now on deriving the corresponding compo- 
nents of the electric field of the ions in the channel. 

Within the channel and at distances, from the ion at 
the origin, smaller than xq , cf. Fig. ^ the electrostatic 
potential is given by Eq. (|10|l . which corresponds, in 
the first order in x/xq to a ID Coulomb potential, cf. 
Eq. (|ll|l . The corresponding components of the electric 
field then formally (expansions of the various functions 
involved have to be taken) read: 



-KQ{R/xQ)jQ{r±/xQ)eTq){-\x\/xo), 



El 



i? < X < xo, r± < R, 



(19) 



El' = - 
i? ^ X <C xo. 



rifo(i?/xo)Ji(r_L/xo)exp(-|x|/xo), 



r I <R. 



(20) 



These expressions display the role of the image force: Jq 
is unity at the origin, so the longitudinal field is constant, 
non-zero and positive. On the other hand, Ji vanishes at 
the origin, so there is no component of transverse force on 
the charge. The ID Coulomb force, which strongly pre- 
vents more than one ion to be inside the pore, therefore 
reads: 



F^eEl^ix 



= 0) = 



.(log(2xo/i?)-7)- (21) 



Another particular sector, of the electrostatic potential, 
corresponds to the cylindrical shell surrounding the pore, 
at i? <C 7"^ ^ Xq, still within the distance xg along the 
channel, from the origin. In this region, Eq. IjKil) shows 
that the potential is 2D-hke. With the help of Eq. l(TF|l . 
which includes higher order corrections, the correspond- 
ing components of the electric field then formally read: 



j^out 



rA:o(r_L/xo)exp(-|x|/xo), 



-^out 



ri4:i(r_L/xo)exp(-|x|/xo), 



(22) 



(23) 



Finally, the last sector corresponds to the region outside 
the cylinder of radius xq and length 2xo where the po- 
tential is 3D-like and isotropic (cf. Eqs. ©, l|14|l and 
(|17|l ). This potential is still unusual, within the water 
channel, because, due to the effect of image charges, it 
is screened by the membrane only: </5(r) = e/emf, ( x > 
Xq, V rj^ ), ( r_L > Xq, V X ). The electric field then 
corresponds to the usual 3D electric field screened by the 
dielectric constant of the membrane. We denote this field 
by: 



^GXt 



( X > xo, V ), ( > Xq, V x ). 



(24) 
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V. IMAGE SELF-ENERGY OF AN ION IN AN 
INFINITE CHANNEL 



The first term in Eq. H30|l . is associated with the electric 
field inside the cylinder, cf. Eqs. H19|) and (|20|1 . and reads: 



The total self-energy of a charge in the channel may 
be derived, with the help of the self-potential, by the 
identity: 



= -^(0). 



(25) 



where the upper script (oo), refers to the fact that an 
infinite channel is considered. As for the electrostatic 
potential, in Eq. H25(l , we consider the self-energy arising 
due to image charges, the image self-energy, and measure 
it with respect to the Born energy of the charge e^/2eu,a, 
where a is the radius of the charge. As for the electro- 
static potential, we decompose this image self-energy, by 
considering a part, Si, arising from the region x <^ xq, 
and a part related to larger distances, Sg""-*, where the 
potential is the regular 3D-like (but screened everywhere 
by the membrane). Hence: 



J](oo) 



^(oo) 



(26) 



We focus first on the part arising from the ID 
Coulomb, Suppose that the charge is brought to 

the origin along the channel. The corresponding self- 
potential is given by Eq. (|10|l at a; = and r± — 0. The 
corresponding self-energy then reads: 



(oo) 



1 



2 e^xo 



(log(2xo/i?) - 7) • 



(27) 



This self-energy may be evaluated by a similar way by 
considering that the charge has been brought to the ori- 
gin from the membrane {i.e. perpendicular to the pore). 
Then, Eq. H18|l has to be used for the self-potential which 
brings us back to Eq. |[77|) . 

In order to determine the total image self-energy, and 
compare it with exact numerical results, we need to take 
care of boundary terms, arising from the cylindrical ge- 
ometry of the system and which influence the final value 
of the self-energy. We find it therefore convenient to re- 
express Eq. (|25() in the form: 



^(oo) ^ e dVp(r)v3(r), p(r) = 5(r), 



(28) 



where V is the total volume of the system. Using the 
inhomogeneous Poisson equation, we may express the 
density of test charge as a function of the potential, in 
Eq. ||^, to obtain: 



1 

8^ 



dVe{r){Vip{r)f 



(29) 



We focus first on contributions to the image self-energy, 
within Xq. 



Stt Ji„ 87r 7o„t 



(30) 



(oo) 



^lAn — 



Ae^XQ xq Jo 



^o(-) / [JoiU) + Jiiu]m 



where = r±/xo. Expanding the integral in R/xq 
leads, in the lowest order in R/xq, to: 



,(oo) _ 1 e^e'^R'^ 



•r^^,oo; _ 

^l.in — 



2 4e2 



K^{R/xo), 



(32) 



where higher order terms in R/xq are neglected. Finally, 
using Eq. 0, yields: 



sS = 7 7^(log(2W^)-7), (33) 



4 emXo 



With the same notations, the second term, in Eq. H3()|l . 
originates from the electric field within xq but outside 
the cylinder, cf. Eqs. (|22|l and H23|) . and reads: 

^tlt = 7^ /' di±U [Km + Km] ■ (34) 

if^mXo Jr/xo 

In the lowest order in R/xq, the result is: 

„2 



(35) 



where the second term is a boundary correction, due to 
the cylindrical geometry. Including this boundary correc- 
tion, the image self-energy, in the lowest order in R/xq, 
reads: 



t(00) 



1 



2 ernXo 



(log(2xo/i?) - 7) - 0.13- , (36) 

ZemXo 



where the first term returns us back to Eq. (|27l) . 

We now proceed to determine the part of the self- 
energy, Sg"""*, at distances x ^ xq. This reads: 



Y^(oo) _ £m 
2.3 - 



Jext 



extN2 



(37) 



where the external volume is beyond the cylinder of 
length 2xo and radius xq. With the help of Eq. I|24|l . 
this reads: 



^(00) 



26j7^Xq 



(38) 



and should be contrasted with the naive expression: 
Sg"""* = e^/2emXo, based on a spherical symmetry. 

From Eqs. H36|) and (|38|l . the total image self-energy, 
Eq. H26() . to the lowest order in R/xq, reads: 

= (log(2:ro/i?) - 7) + 0.26-^. (39) 
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VI. NUMERICAL ESTIMATES FOR THE 
INFINITE CHANNEL 

The most straightforward and accurate way, of deter- 
mining the value of the self-energy of an ion, Eq. (|25|l . 
has been taken by Parsegian who has used a formula™, 
similar to Eqs. and which has been substituted 
in Eq. H25I) . and evaluated numerically. With the help of 
our Eqs. (P) and |(2Jl, Eq. if^ reads: 



TABLE I; Image self-energy, of an ion in an infinite chan- 
nel, in various units from the exact numerical calculation, see 
Eqs. (0, and discussion below them as well as Ref. Q, 
and from our asymptotic formulas, Eq. (I39II . We have taken 
em ~ 2 and e^, = 80. 



J](oo) 



dk 



)Ko{kR)Ki{kR) 



,Ko{kR)h{kR) + eraKiikR)IoikR) 



Self-energy unit 


Exact 


Eq. 




e'/tmR 


1.08 


1.11 




0.17 


0.17 






6.80 


7.01 






16.0 


16.6 


= 3.0A) 



where the energy is measured with respect to the Born 
energy, and may be evaluated with any standard soft- 
ware. We shall refer to this solution, for the image self- 
energy, as the exact numerical solution. 

Our analysis of the previous sections, on the other 
hand, has revealed some important features of the elec- 
trostatics of ion channels, in particular the crossover 
length xq, the ID nature of the Coulomb potential and 
expressions for all quantities of interest, including the 
self-energy of Eq. (|4U|I and the main contribution to 
this image self-energy, the ID Coulomb part. We pro- 
ceed to evaluate these expressions, by taking the case 
of, e.g. the model-pore gramicidin A which has a radius 
R = 3.0A and a length L = 25. OA, and for the typical 
room-temperature values: = 2 and £„, = 80. 

The crossover length, Xq, has to be determined self- 
consistently from Eq. |7J • In the logarithmic approxima- 
tion, this length reads: 



xo « i?W-^V0.51og(2e^/e„0 - 7. (41) 

V ^^m 

Using a more precise graphical method, cf. Fig.|21 yields: 



Xo ~ 6.3i?, 



(42) 



where we have used the fact that e™ = 2 and e„, — 80. 
In the case of the model-pore gramicidin A {R — 3.0A): 
Xq « 19. OA. The peculiar potential of Eq. Hll|l plays a 
dominant role over this length, from the origin. Beyond 
the length xq, the ID Coulomb potential crosses over to 
the 3D one. We have therefore to compare the length of 
the channel, L = 25. OA, to the region where the anti- 
confinement dominates: 2xo « 38. OA. 

Obviously, such a value of 2xo > L shows that the ID 
Coulomb potential dominates in the model nano-pore. 
We shall come back on this crucial fact, in the next sec- 
tion, when dealing with finite-size effects. Substituting 
the expression of xq, in Eq. H10|) . leads to: 



<f{x) « 



v/2 (log(2xo/i?) - 7) 



exp(. 



\x\/xo), 
(43) 



which displays the "geometric mean" screening: oc 
yj^w^mi and is affected by a logarithmic factor charac- 
teristic of the intrinsic ID nature of the system (see also 



Ref. where the geometric-like screening has been con- 
sidered qualitatively). 

As a consequence of Eq. (|43|) , if an ion is in the pore, 
any additional ion will feel a constant repulsive force, cf. 
Eq. Substituting the expression of Xo, cf. Eq. Q, 

in Eq. (|^ . this force simply reduces to: 



F 



2e2 



(44) 



This is an unrealistically large (because we consider an 
infinite channel) force: Fa ~ 48.9 eV/A. The strong re- 
pulsion originating from Eq. H44() is therefore responsible 
for the large image self-energy of the charge, cf. Eq. H27|l . 
which, with the help of Eq. {Tj), and to the lowest order 
in R/xq, reads: 



^(oo) 



2 \/ f^w^rnR 



V2(log(2xo/i?)-7)-Cg, (45) 



where Cg = 0.13(e^/2emXo), is the geometric correction, 
appearing in Eq. (|36|) . With the help of = 2 and 
= 80, Eq. igHl) yields: 



^(00) 
where C, 



(44.02/i?- 2.90/i?)fcsT w (41.12/i?)/ci3T, (46) 



g ~ {2.90/ R)kBT. This correction is not negli- 
gible when it comes to compare our results with precise 
numerical estimates. For R = 3A, Cg ~ 0.97fcsT, and 

this yields: E 



(00) 



crease results: 



IS.lksT. For i? = 2A, a drastic in- 
« 20.56A:bT, where Cg « LAbknT. 
On the other hand, for R — 5A, this quantity is reduced 
to: « 8.22fcsT, Cg « 0.58fcsT. 

To estimate the total image self-energy one should con- 
sider Eq. (|39|l . which takes into account of the large- 
distance part of the energy. Adding this contribution, to 
Eq. igni), yields: 



S(°°) w {49.82/ R)kBT. 



(47) 



For R ^ 3A, Eq. yields: « 16.6A:bT. Numeri- 

cal estimates for in various units, are displayed in 

Tab. m They agree with the exact results of numerical 
simulations. 

Notice that, to express the self-energy in units of 
ksT, where fc^ is Boltzmann's constant, we have found 
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FIG. 6: Schematic view of a long channel, for which: L 3> 
2x0, of. Eq. Q. 



it convenient to introduce Bjerrum's length, Is- At 
room temperature and in water, Bjerrum's length reads: 
Ib = e"^ /ewksT. To be able to compare safely our results 
with previous studies^ we take the value Ib ~ 7.038A 
(Jordan converts the values of Levitt (in ksT) in units of 
e/e^R and this value of Ib agrees with the conversion). 



VII. FINITE-SIZE EFFECTS: 
TRANSLOCATION ENERGY IN 
NANO-CHANNELS 

In the case of a finite-length channel, calculations sim- 
ilar to the infinite length case lead to coupled non-linear 
integral equations, at the level of the boundary conditions 
imposed on the electrostatic potential of the inhomoge- 
neous Poisson equation. Even an approximate solution of 
these boundary equations seems intractable, analytically, 
and numerical simulations appear to be necessary. Such 
a numerical task has been undertaken in Refs. and '3, 
by two different numerical methods, showing that the 
self-energy becomes of the order of E'^^-' « 6.7kBT 

for a length L = 25. OA and a radius R = SA. The latter 
is defined as: 

= (e/2)(^(0)-^(L/2)), (48) 

which is simply the work required to bring the charge 
from the entrance of the channel, at a; = L/2, to its 
center, at a; = 0. Two formulas were then proposed to 
take into account of the finite size effects. 
For large lengths Parsegian has proposed^: 

E(^)«S(-)-^log(^5-), (49) 

where the first term is the self-energy of the infinite chan- 
nel and the second term corresponds to the self-energy of 
a charge in a membrane without a pore. On the light of 
the previous sections, this "large length" ansatz, should 
work for lengths larger than 2xa. We may understand 
this result on the basis of Fig. IHl by decomposing the 
self-energy into two parts, while bringing the charge from 



FIG. 7: Schematic view of a nano-length channel, for which: 
L <^ 2xo, cf. Eq. O. This case corresponds to the typical 
channels we are interested in, e.g. gramicidin A for which 
L ^ 2.5nm, i? « SA and 2xo « 12.6i? ^ 3.8nm. 

infinity to the origin: 

sW^si-) + ^(l-|), (50) 
2e„ xq L 

where the first term corresponds to the ID Coulomb part 
of the self-energy of the infinite channel and the second, 
to the energy required to displace the charge from L/2 
to Xq (for the finite channel image charges do not play 
any role beyond L/2 from the origin and the correspond- 
ing self-energy contribution has been neglected). Up to 
numerical coefficients, the first two terms yield and 
the last term yields the (reducing) effect of the mem- 
brane. Unfortunately, Eq. H49() seems to work only for 
unrealistically long channels, well beyond 2xo, i.e. the 
percentage of error, with respect to numerical simula- 
tions, becomes less than 5% for lengths larger than 40R. 

For the most relevant case of small length channels, cf. 
Fig. [7| Jordan has proposed, without any justification 
other than curve-fitting, the following expression^: 

E(^) « E(°°)[l-exp(-6(5)], (51) 

where b — 0.131 and S is the reduced channel length, 
S = L/2R . We see immediately that the characteris- 
tic length appearing in the exponential corresponds to: 
Lq ~ 2.44x0. The "small-length" case will therefore obvi- 
ously correspond to the regime where L <ti 2xo. We now 
proceed to justify the curve-fitting Eq. (|51|1 and improve 
on the result of Jordan. 

As is clear from the previous lines, the knowledge of 
Xq allows us to consider that channels of characteris- 
tic length: L ^ 2xo, are finite, with respect to the 
ID Coulomb part of the self-energy, whereas channels 
of characteristic lengths: L 3> 2xo, are formally infinite 
with respect to the same quantity. As a matter of fact, 
the ID Coulomb part of the self-energy of a channel of 
length L ^ 2xq saturates: sj^"* « On the other 

hand, this self-energy, for a channel of length L <C 2xo, is 
determined with the help of Eq. igSl and Eq. gHl) • The 
first term in Eq. H48|) corresponds to and the second 
gives the leading correction in L/xq. Actually, higher or- 
der corrections in L/2xq are important. As in Eq. 143|l . 
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such corrections may be re-summed in an exponential 
factor. This yields: 

«I](°°^[l-exp(-i/2a;o)], L<2a;o, (52) 

where S^""-* is given by Eq. H45|l . which justifies the func- 
tional form of Jordan's expression on a physical basis. In 
the lowest order in L/2a;o, '^[^\ in Eq. H52|l . increases lin- 
early with L, a feature of the ID nature of the Coulomb 
potential: T,[^^ « E^°°^L/2xo. For the case of interest 
{em = 2 and = 80), the image self-energy of a channel 
of length L and radius R reads: 

/r\ 44 02 

R) « (-^ - Cg)[l - exp[-V6.3]] fc^T, 

L < 2x0, (53) 

where S = L/2R and Cg w {2.90/R)kBT from the previ- 
ous section. Actually, Eq. ifS^ works better than Eq. 
at lengths L ^ 2xo. This implies that, in Eq. (|49|l . cor- 
rections due to the presence of the membrane beyond 
\x\ > xo manifest only for very large lengths, L ^ 2xo, 
which is beyond the lengths of the channels we are inter- 
ested in. 

We also notice from Eq. H53|l , that the self-energy can- 
not be expressed in terms of a single scaling parameter, 
here 6 = L/2R. Eq. H53|) depends separately on L and 
R. This fact is very clear from the work of Levitt, that 
we have reproduced in Tab. ^] For a given parameter S, 
e.g. 6 = 6.25, the channel (L = 25A, R = 2A) is char- 
acterized by E^^' « 13. BfcsT whereas the channel (L = 
37.5A, R = SA) is characterized by: E^^) « S.gifcsT. 
This is a drastic energy variation, which is well repro- 
duced by Eq. 

Eq. I|53|) gives therefore a fit-parameter free expression 
for the self-energy of a charge in nano-channels, in the 
range of validity: R <^ L <^ 2xq. As can be seen from 
Tab. m the average percentage of error (in the regime 
of validity) of Eq. (|53|) . with respect to numerical sim- 
ulations, is: 2.16% for the data of column DL and PJl 
and of 4.14% for column PJ2. With respect to the data 
of Levitt^ we are within 3.% of error and with respect to 
the data of Jordan^ within 5.%, i.e. similar than the er- 
ror bar of Jordan's curve-fitting formula. Notice also that 
the accuracy extends beyond L = 2xq (or 5 = 6q = 6.3), 
with an error still within 3.%, with respect to PJl and 
within 6.%, with respect to PJ2. In the limit L 3> 2xq, 
Eq. (|53|l also gives a better accuracy than Eq. H49|) the 
values of which have been displayed in the last column 
of Tab. ^ Better accuracy, with respect to the numeri- 
cal simulations, may be reached if Cg is considered as a 
fitting parameter but we shall not dwell on this point as 
our fit-free results have already a very good accuracy. 

VIII. CONCLUSION 

In summary, we have shown that finite-T transport, 
across ionic channels, is determined by zero-temperature 



(as there is no screening^, by counter-ions within the chan- 
nels we are interested in) classical electrostatics. Our 
main results have shown that ions in the water chan- 
nel of a lipid membrane interact via a one-dimensional 
Coulomb potential, cf. Eq. Hll|l . below a crossover length 
xq, cf. Eq. Q, and a three-dimensional Coulomb poten- 
tial beyond this length. The length xq has been shown 
to be of the order of the length of, e.g. the model pore 
gramicidin A. An expression for the image self-energy of a 
finite nano-channel has then been obtained, cf. Eqs. 1)52(1 
and l|45|l . with no fitting parameter and an agreement 
within an error of 5% with the existing numerical simu- 
lations available from the literature on the subject. Com- 
bining Eqs. H52|) and H45|) the most general expression for 
the image self-energy of a channel of length L and radius 
R reads: 

E(^)(i,i?) = 

^ ^^p V2(log(2xo/i?)-7)[l - exp(-L/2xo)], 

R<^L<. 2a;o, (54) 

where 7 « 0.577 is Euler's constant, x^ is determined 
by Eq. Q and A/" = 1 — 0.13, is a numerical coefficient 
originating from Cg, cf. Eq. (|45|) . The latter may be of 
the order of a few fc^T (depending on the characteristics 
of the channel under consideration) and has to be taken 
into account for numerical estimates. In the lowest order 
in L/2xq, Eq. H54|) increases linearly with the length of 
the channel, a feature of the ID Coulomb potential. 

The image self-energy of Eq. ((54|l corresponds to 
the translocation, or activation, energy through a finite 
length channel, e.g. the energy required to displace a 
charge from the entrance of a channel to its center. In 
relation with experiments, the existence of such elec- 
trostatic barrier requires a threshold voltage difference, 
eAT4(L,i?) = Y.[^\l,R) (where the self-energy is in 
units of eV), across the channel in order to pass an ion 
through it (no chemical potential difference is taken into 
account here). Equivalently, this corresponds to the ex- 
istence of a threshold electric field: Ec = AVc{L,H)/L, 
applied along the length of the channel. As an inverse 
problem, measurement of this threshold field would then 
give access to information on the channel, its radius or 
length, with the help of Eq. 

The image self-energy of Eq. ((54|l determines the main 
exponential dependence of the conductance of a nano- 
scale channel: 

G (xexp[-J:[^\L,R)/kBT]. (55) 

Up to now, all our estimations where performed at 
room temperatures for — 80. Actually, the di- 
electric constant of water is sensitive to temperature 
and varies from e„(273 K) = 87.78 to e^,(323 K) = 
69.91. Basic estimations of Ei = Ei(T) which take 
into account of this temperature variation of e„ show 
that the ratio Si (T)/fcBT becomes approximately tem- 
perature - independent (the coefficient Ct defined as: 
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TABLE II: Image self-energy of an ion in a channel of length L, in units of fcsT, as a function of the reduced channel length 
5 = L/2R, where R is the radius of the channel. We have taken = 2 and = 80. Energies are in units of ksT and 
lengths are in units of A. The crossover length, Lq = 2xo between large and small channels, corresponds to the universal 
value (independent of the length and radius of the channel): Sq — 6.3. Our asymptotic calculations are valid for lengths: 
R < L < 2x0. Eq. H53|l . is taken with no fitting parameter. The percentage of error appearing in columns DL, PJl and PJ2, is 
that of Eq. I53II with respect to the numerical simulations. 
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"As determined from Table I of Levitti. 

'As determined from Table I, column 2, of Jordan™. 

■^As determined from Table I, column 3, of Jordan™. 



^i{T)/T\t=27zk = Ct^i{T)/T\t=3,2zk is equal to 
Ct ~ 1.18 for a temperature-independent dielectric con- 
stant and to Ct ~ 1.00 vifith the above temperature- 
dependence of Eu)). This implies that, despite the fact 
that the conductance is affected by an exponentially 
small factor, its thermal dependence is weak so that the 
ionic flow through channels of cell membranes is sta- 
ble, being approximately constant over more than two 
decades of units of temperature, around room tempera- 
ture. 
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